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One-band Hubbard model with transverse anisotropy is considered at density of electrons n = 0.4. 
It is shown that when the anisotropy is appropriately chosen, the ground state is ferromagnetic with 
magnetic order perpendicular to the anisotropy. The increasing of the ratio y , where t is the hopping 
parameter and U is the Coulomb repulsion, decreases the Curie temperature, and the system arrives 
at the quantum critical point (Tc = 0). The result is obtained introducing Schwinger bosons and 
slave Fermions representation of the electron operators. Integrating out the spin-singlet Fermi fields 
an effective Heisenberg model with ferromagnetic exchange constant is obtained for vectors which 
identifies the local orientation of the spin of the itinerant electrons. The amplitude of the spin 
vectors is an effective spin of the itinerant electrons accounting for the fact that some sites, in the 
ground state, are doubly occupied or empty. Owing to the anisotropy, the magnon fluctuations drive 
the system to quantum criticality and when the effective spin is critically small these fluctuations 
suppress the magnetic order. 

PACS numbers: 64.70.Tg, 71.10.Fd, 75.50.Ee, 75.40.Cx 



Quantum phase transitions (QPT) arise in many-body 
systems because of competing interactions that sup- 
port different ground states. At quantum critical point 
(QCP)the matter undergoes a transition from one phase 
to another at zero temperature. A nonthermal control 
parameter, such as pressure, doping or magnetic field, 
drives the system to QCP. Quantum phase transitions are 
a subject of great interest jlJ-Q . At this point, the quan- 
tum critical fluctuations give rise to unconventional tem- 
perature dependence of magnetic, thermal, and transport 
parameters [JQ. 

Quantum phase transition can be induced in a wide 
range of materials. Most prominent experimental realiza- 
tion of ferromagnetic quantum phase transition is found 
in magnetic properties of LiHoF^ Q. At low temper- 
ature the magnetic degrees of freedom in this material 
are the spins of the holmium atoms. They have an easy 
axis, and below Tq = 1.53K the compound is ferromag- 
net. In Q the authors measured the magnetic order as 
a function of temperature and a magnetic field applied 
perpendicular to the easy axis. The increasing of the 
magnetic field reduces Curie temperature monotonically. 
When it is larger than some critical field (about 50fcOe) 
the long-range order in the material is destroyed even at 
zero temperature. The spin flip operators are the quan- 
tum fluctuations which drive the system to the quantum 
critical point 0, ■ Another examples of field-induced 
quantum phase transitions are discussed in the review 
articles @,HH. 

Pressure driven transformation of GeRu^Gei from fer- 
romagnet into paramagnet, at zero temperature, is stud- 
ied in 0] . The suppression of magnetic order at the crit- 
ical pressure p cr = 67kbar is accompanied by non-Fermi- 
liquid behavior. 

The magneticproperties of U Rh\- x Ru x Ge alloys are 
investigated in (9|. The Curie temperature vanishes lin- 
early with "x" and the ordered moment is suppressed in 



a continuous way at x cr = 0.38. The thermal, transport, 
and magnetic properties of U Rh\^ x Ru x Ge near the crit- 
ical concentration are investigated. The data provide ev- 
idence for continuous ferromagnetic quantum phase tran- 
sition. 

The earliest theory of a ferromagnetic quantum phase 
transition was the Stoner theorv[lG|. In later paper [l[ 
Hertz derived an effective Ginzburg-Landau- Wilson the- 
ory starting from a microscopic theory of itinerant elec- 
trons with four-fermion interaction. He concluded that 
in the dimensions d — 2 and d = 3 the quantum phase 
transition from paramagnet to fcrromagnet in isotropic 
system is a second-order transition with mean-field crit- 
ical behavior. In contrast to this it was shown 
that the quantum phase transition from a metallic para- 
magnet to an itinerant ferromagnet in d = 2 and d = 3 is 
discontinuous. It was proved that this conclusion is true 
and for anisotropic systems with magnetic order parallel 



to anisotropy axis [14 1 



In the present paper the three-dimensional Hubbard 
model with transverse anisotropy is investigated. The 
experimental observations in LiHoF± inspire that the 
anisotropy could be of grate importance for the exis- 
tence of a quantum critical behavior in the system. It 
is shown that when the anisotropy is appropriately cho- 
sen, the ground state is ferromagnetic with magnetic or- 
der perpendicular to the anisotropy. The increasing of 
the ratio jj, where t is the hopping parameter and U 
is the Coulomb repulsion, decreases the Curie temper- 
ature. Owing to the anisotropy, the spin flip fluctua- 
tions (magnons) drive the system to quantum criticality 
and when the effective spin, which accounts for the fact 
that some sites, in the ground state, are doubly occupied 
or empty, is critically small these fluctuations suppress 
the magnetic order. Critically high double occupancy is 
phenomena of basic relevance to second-order quantum 
phase transition in itinerant magnets. 
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We consider a theory with Hamiltonian 

h = - t E (c+ Cj CT + /i.e.) + U E - /^E n, 

- ./')^'.sf.s; • /r.sf.sf • s;s]' (1) 

where c^. and Cj CT (cr =f, J,) are creation and annihilation 
operators for spin-1/2 Fermi operators of itinerant elec- 
trons, rii a = c^Cia, rii = rii\ + n±, t > is the hopping 
parameter, U > is the the Coulomb repulsion, and /i 
is the chemical potential. The exchange constant in the 
Heisenberg term is ferromagnetic J' > 0, k is a param- 
eter of anisotropy, and the spin of the itinerant electron 
is 

s^lJ2^ c ^ (2) 

aa' 

where t x ,t v ,t z are the Pauli matrices. The sums in 
Eq. ([1]) are over all sites of a three-dimensional cubic lat- 
tice, and denotes the sum over the nearest neigh- 
bors. 

We represent the Fermi operators, the spin of the itin- 
erant electrons and the density operators in terms 
of the Schwingcr bosons ((pi, a ,(p^ " ) and slave Fcrmions 

(hi,h^~ ,di,df). The Bose fields are doublets (a = 1,2) 
without charge, while Fermions are spinless with charges 
1 (di) and -1 (hi). 

c lt = h+cpn + (pf 2 di, Cii = hfip i2 -tp^di, 

Ui = l- Jlfhi + d+d, , ^ = -^2 <fit T aa> Via' , 

aa' 

4\ c ^ c tx c ^ = d i d * ( 3 ) 

vtiVii + vt 2 H>i2 + dfdi + hfhi = 1 (4) 

To solve the constraint (Eqj4j , which entangles the 
Fermi and Bose operators, one makes a change of vari- 
ables, introducing Bose doublets Q a ((+„) [H| 

i 

fia = Qa (l - h~l hi - dfdi) 2 , 

<p& = ai-hihi-dfd^ . (5) 

Now, only the new Bose fields are constrained Clada = 1- 
In terms of the new fields the spin vectors of the itinerant 
electrons have the form 

S? = \ E COC^ [1 - h+hi - d+di] (6) 

aa' 

When, in the ground state, the lattice site is empty, the 
operator identity hi hi = 1 is true. When the lattice 
site is doubly occupied, dldi = 1. Hence, when the lat- 
tice site is empty or doubly occupied the spin on this 



site is zero. When the lattice site is neither empty nor 
doubly occupied [hlh = dldi = 0), the spin equals 
si = l/2ni, , where the unit vector 

<=EC^'(»' (n? = l) (7) 

aa' 

identifies the local orientation of the spin of the itinerant 
electron. 

The Hamiltonian Eq.([T]), rewritten in terms of Bose 
fields Eq.© and slave Fermions, adopts the form 

(ij) 

+ (d+h+ - d+h+) (CnOa - CaGi) + he] 
x (1 - hlh - did,) ~ 2 (1 - hfhj - d+dj) ~ 2 
+ U E 4di - fj, E (1 - h+h + did,) 

i i 

- j'Yl\. n >j+ Kn X+ n >j] ( 8 ) 

m 

x (l - hlh - dfdi) (l - hjhj - d+dj) . 

An important advantage of working with Schwinger 
bosons and slave Fermions is the fact that Hubbard term 
is in a diagonal form. The fcrmion-fcrmion and fcrmion- 
boson interactions are included in the hopping and spin 
exchange Heisenberg terms. To proceed we approxi- 
mate the hopping term of the Hamiltonian Eq.® setting 

(l — hlh — dldlj 2 ~ 1 and keep only the quadratic, 
with respect to Fermions, terms. This means that the 
averaging in the subspace of the Fermions is performed 
in one fermion-loop approximation. Further, we repre- 
sent the resulting Hamiltonian as a sum of two terms 

h = h + h int , (9) 

where 

ho= - t E {dfck - h+h + he.) d t d ^ 

(ij) 

- n^2(l-hthi + dtdi), (10) 

i 

is the Hamiltonian of the free d and h fermions, and 

h int = - <E Mdi - h+h) (ClaQa ~ 1) (11) 

(ij) 

+ ( d ihl - d+h+) (G1O2 - Q2Q1) + he] 

- J'E[< n i + Kn K +n Hl (12) 
(ij) 

x (1 - hlh - did,) (1 - hjhj - d+dj) 

is the Hamiltonian of boson-fermion interaction. 

The ground state of the system, without accounting for 
the spin fluctuations, is determined by the free-fermion 
Hamiltonian ho and is labeled by the density of electrons 

n = 1- < h+hi > + < dldi > (13) 
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(see equation ([3]) ) and the " effective spin" of the electron 
1 



m 



(1- < hfhi >-< d+d, >) . (14) 



At density of electrons n = 0.4 one calculates the chem- 
ical potential /i as a function of t/U and utilizes the result 
to calculate the effective spin of the itinerant electron 
"m" as a function of the control parameter t/U . The 
result is depicted in Fig. ([1} -black triangles-left scale. 

Let us introduce the vector, 



M? 



Ml 



171 . 



(15) 



It is convenient to rewrite the effective hamiltonian in 
the form 

h eff = ~jJ2l M i M ! + rMfMj + MfMj], (20) 

m 

with exchange constant J = J" + J' and effective param- 
eter of anisotropy r = (J" + kJ')/{J" + J"). 

The exchange constant J, the effective spin m and the 
effective anisotropy parameter r are functions of the ratio 
t/U. They are calculated at density of itinerant electrons 
n = 0.4, for anisotropy parameter n = —0.83 and J' jt = 
10. The functions are depicted in Fig. ([I]). 



Then, the spin- vector of itinerant electrons Eq. © can be 
written in the form 

S, = ^-Mi (1 - h+ hi - dj di) , (16) 

where the vector identifies the local orientation of 
the spin of the itinerant electrons. The contribution of 
itinerant electrons to the total magnetization is < Sf >. 
Accounting for the definition of m (see EqfT3]), one ob- 
tains < Sf >=< Mf >. 

The Hamiltonian is quadratic with respect to the 
Fermions di , df and hi , hf , and one can average in the 
subspace of these Fermions (to integrate them out in the 
path integral approach). As a result, one obtains an ef- 
fective model for vectors with Hamiltonian 

m 
(«) 

The effective exchange constant J" is calculated in the 
one fermion-loop approximation, in the limit when the 
frequency and the wave vector are small. At zero tem- 
perature, one obtains 



J = 



t 1 

6m 2 N 



2t 2 



(J2cosk u )[0(-ei)-e(-4)] (18) 



3m 2 U N 



1 - 9(-e! 



o(-4)] 



where N is the number of lattice's sites, e\ and e'/. are 
Fermions' dispersions, 

e\ = 2t(cos k x + cos k y + cos k z ) + fj, (19) 
= — 2t(cos k x + cos k y + cos k z ) + U — fi, 

and the wave vector k runs over the first Brillouin zone 
of a cubic lattice. The "m" dependence of J" is a con- 
sequence of the definition of the vector M, Eq. (fT5"|) . and 
J"m 2 doesn't depend on "m". 
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FIG. 1: (Colour on-line) The effective spin of the itinerant 
electrons m as a function of the control parameter f/[/-black 
triangles(left scale) The exchange constant J/t as a func- 
tion of t/U-red rhombuses(right scale). Inset: The effective 
anisotropy constant r as a function of t/U. The functions 
are calculated at density of itinerant electrons n = 0.4, for 
anisotropy parameter k = —0.83 and J' /t = 10. 

For chosen parameters the exchange constant is posi- 
tive J > and a ferromagnetic order along z axis, perpen- 
dicular to the anisotropy, is stable. To study this phase 
one introduces the Holstein-Primakoff representation of 



the spin vectors Mj(a^~, a 3 ) 



Mf = Mj + /.I/; = ^2m - a+a 3 aj 



Mr = Mf - iMV = af Jim - a+a 3 



M 



m — a ■ a 7 -, 
j J 



(21) 



where af , aj are Bose fields and m is the effective spin of 
itinerant electrons. In spin-wave approximation, in mo- 
mentum space representation the effective Hamiltonian 
fEq |2"0)) adopts the form 



h 



eff 



E [ £ k a t a k ~ 7fc (flfe at k + a fe a_ fc )] 



4 



7 fc = Jm 



Jm [6 — (1 + r) (cos k x + cos k y + cos k z )\ 



(cos fc^ + COS fey + COS fc z ) . (22) 



To diagonalize the Hamiltonian one introduces new Bose 
field afc, a^: = Ufe ctk + «fc ct^,, where the coefficients 
Uk and Vk are real functions of the wave vector. The 
transformed Hamiltonian adopts the form 



(23) 



with dispersion 

E k 
h 



2Jm ^(3 - 4) (3 - rS k ) 
cos fc^ + cos k y + cos fc z 



(24) 



and vacuum energy E® = i [E^ — 

The dispersion Eq. (|24)) shows that the ferromagnetism 
with magnetic order perpendicular to the anisotropy axis 
is stable if the effective anisotropy constant satisfies the 
condition — 1 < r < 1. The inset in Fig.([T]) shows that 
all calculated values of the effective anisotropy constant 
satisfy the above condition. 

The dispersion is equal to zero at k = (0, 0, 0). There- 
for, a^-Boson describes long-range excitation (magnon) 
in the system. Near the zero vector the dispersion adopts 
the form Ek oc c s |k| with spin- wave velocity 



Jmi/ 6(1 — r). 



(25) 



The unusual, for ferromagnetism, dispersion is in con- 
sequence of anisotropy. When the isotropy is restored 
(r = 1) one obtains the well known ferromagnetic disper- 
sion. 

The magnetization is defined by the equation (|21l) . At 
Curie temperature Tc, the magnetization is zero and one 
calculates the critical temperature solving the equation 



2m + 1 



1 \ -» £k 
v 2-^~E k 

k 



N 



1 



(26) 



When the parameter of anisotropy r is fixed, the critical 
temperature decreases, decreasing the effective spin m. 
At quantum critical point (Tc = 0) one obtains a depen- 
dence of the critical value of the effective spin m cr on 
the anisotropy parameter. The relation is given by the 
equation 



2m c 



1 



N 



2N 



6-(l + r)4 
y/{3-S k )(3-rS k ) 



(27) 



with 8k from Eq. (|2"4")) . and is depicted in Fig©. 

The figure (2) demonstrates very well the nature of 
the quantum criticality in itinerant ferromagnets with 
anisotropy. At quantum critical point the effective spin, 
which is the fermion contribution to the magnetization 



without accounting for the spin fluctuations, is nonzero. 
Magnon fluctuations suppress this magnetization to zero. 
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FIG. 2: (Colour on-line) The effective spin m cr , at quan- 
tum critical point, as a function of the effective parameter 
of anisotropy. 



Hence, the magnon (spin flip) fluctuations drive the sys- 
tem to quantum critical point. This is a second order 
quantum phase transition from ferromagnetic phase with 
long-range magnon fluctuations to paramagnetic phase 
with gapped magnons. 

This is in contrast to quantum phase transition in 
isotropic itinerant ferromagnets (r = 1). For this sys- 
tems the quantum phase transition is from ferromagnetic 
phase with long-range magnon fluctuations to param- 
agnetic phase without spin- flip fluctuations (m cr = 0). 
This explains the non-second order nature of the transi- 
tion. 

Utilizing the dependence of the effective spin m and 
the exchange constant J/t on the control parameter t/U, 
see Fig. ([IJ, one can obtain the dependence of the di- 
mensionless temperature Tc/t on the ration t/U. The 
phase diagram in the plane of temperature Tc /t and con- 
trol parameter t/U is depicted in figure ([3]). At density 
of electrons n = 0.4, anisotropy parameter k = —0.83 
and J' /t = 10 the quantum critical value of the ratio is 
t/U = 2. 

In the present paper the Hubbard model of itinerant 
ferromagnetism with transverse anisotropy was studied. 
Increasing the ratio t/U decreases the effective spin of 
itinerant electron which in turn decreases the Curie tem- 
perature. The equations (| 13(1 and (fT4"|) show that at fix 
density of electrons n = 0.4 the decreasing of the effec- 
tive spin increases the density of doubly occupied states. 
The result shows that the quantum critical ground state 
state is a state with high concentration of doubly occu- 
pied states. 

This work was partly supported by a Grant-in-Aid 
DO02-264/18. 12.08 from NSF-Bulgaria. 
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FIG. 3: (Colour on-line) The phase diagram in the plane of 
temperature Tc/t and control parameter t/U. At density 
of electrons n = 0.4, anisotropy parameter k = —0.83 and 
J' ft = 10 the quantum critical value of the ratio is t/U = 2. 
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